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Abstract 

We consider general relativistic magnetohydrodynamics from a charged multifluids point-of- 
view, taking a variational approach as our starting point. We develop the case of two charged 
components in detail, accounting for a phenomenological resistivity, providing specific examples 
for pair plasmas and proton-electron systems. We discuss both cold, low velocity, plasmas and 
hot systems where we account for a dynamical entropy component. The results for the cold case 
(which accord with recent work in the literature) provide a complete model for resistive relativistic 
magnetohydrodynamics, clarifying the assumptions that lead to various models that have been 
used in astrophysical applications. The analysis of the hot case is (as far as we are aware) novel, 
accounting for the relaxation times that are required to ensure causality and demonstrating the 
explicit coupling between fluxes of heat and charge. 
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I. INTRODUCTION 



Magnetic fields are ubiquitous in the Universe, affecting physics across a vast range of 
scales. The relevance of electromagnetism for our everyday experience is obvious. Elec- 
tromagnetic effects are also central to many processes in astrophysics and cosmology. The 
strongest known magnetic fields (above 10^^ G) are found in a subclass of neutron stars aptly 
referred to as magnetars , systems that also form the largest (and hottest!) known su- 



pe.o..ucto. BB. Ma,„e«c fields e,ual, .>eva.. on tKe vast, >a^e. .ale of ent.e 
galaxies, and are likely to have played a role in the early Universe as well |5|-l7|. Understand- 
ing the origin and evolution of electromagnetic fields in their many different guises remains a 
fundamental question for modern science. The literature on the subject is (understandably) 
vast yet some problems remain relatively unexplored. This paper concerns one such 
problem. 

Our aim is to develop a model for resistivity in general relativistic magnetohydrodynam- 
ics. By necessity this forces us to consider a charged multi-fluid system (we obviously need 
charged components in relative motion in order to have a charge current!). This part of the 
problem is quite straightforward; we make progress by marrying the standard variational 



model for electromagnetism 



to the charged fluid version of Brandon Carter's convective 



variational description of relativistic fluids |9 



lol- 



Adding a phenomenological resistivity to 



;he mix is not difficult, either. Combining these ingredients we follow the text-book strategy 
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12| and derive the simplified equations of magnetohydrodynamics. As long as we limit 



the analysis to low velocities (cold plasmas) the results follow readily. We demonstrate this 
for the particular problem of a two-component system, composed either of protons and elec- 
trons or a pa.ir-plasma with positrons and electrons, and compare our results to the recent 



literature 



13 



15|. 



The complexity of the problem increases significantly if we turn our attention to high 
velocities and hot plasmas. One reason for this is obvious: In order to describe a hot 
system we need to allow for the presence of heat fiow. However, the problem of heat in 
relativistic systems is known to be difficult, as a naive implementation inevitably leads to 
causality violation and unwanted instabilities 



16 



19| . We avoid falling into this trap by 



building on a recent model that treats the entropy as an additional "fiuid", which couples 
to the substantial matter components through entrainment 2CH22| . This effect represents 
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the inertia of heat, and leads to the thermal relaxation that is required to ensure causality 
and stability. The presence of this coupling makes the analysis less straightforward, and the 
final results are (obviously) less transparent than in the low-velocity case. However, they 
are also more "interesting" . The more complicated setting allows for a number of additional 
features, most notably a coupling between the heat flow and the charge current. From a 



fundamental point of view one would expect such a thermo-electric coupling [23|, but this 
effect has nevertheless not been previously discussed in a relativistic context. 

II. CHARGED RELATIVISTIC MULTI-FLUID SYSTEMS 

This section sets the stage for the discussion by bringing together and adapting established 
results from the literature. The key building blocks are obvious: We need a framework for 
discussing electromagnetism in general relativity, and in order to understand the nature of 
the associated current we also need a multi-fluid formulation for charged components. The 
first part can be found in many standard text-books (see for example 8|). The multi- fluid 
part is less mainstream fare, but the required formalism (mainly designed by Carter and 



colleagues, see jo]. 



lo| for reviews) has been developed to the required level. The marriage 



of the two systems has not been discussed extensively in the literature but, as we will see, 
it is comfortable. 



A. Variational multi-fluid dynamics 

Multi-fluid dynamics arise whenever a system has several components, each in the "fluid 
regime" , which retain their identity. The archetypal such system, known to be well described 
by a two- fluid model, is superfluid ^He 2J]. In principle, one can imagine systems where 



the mean-free path due to inter-species scattering is much larger than that for intra-species 



scattering 25|]. On intermediate scales one can then meaningfully discuss different fluid 
components. This set-up may seem somewhat artificial, but there clearly are systems in 
nature where this separation of scales occurs. One reason why superfluid systems tend 
to require a multi-fluid approach is that the relevant scale deciding the "size of the fluid 
elements" is not the mean-free path (since particle scattering is suppressed in a superfluid) 
but the coherence length of the relevant condensate. This scale is usually much smaller 
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than the mean-free path in the corresponding system at temperatures above the superfluid 
transition, so the system ends up acting as a fluid on much smaller length scales than usual. 
In an astrophysical context, the modelling of mature neutron star cores must account for 
superfluidity (and superconductivity!). Indeed, most applications of the general relativistic 



multi- fluid formalism have been in that problem area, see for example [26 1. 

The model we consider builds on the convective variational principle developed by Carter 
9|. This method deals, in a natural way, with the fact that a variational derivation of the 
equations of fluid dynamics must be constrained. The development takes as starting point 



a Lagrangian for the matter. A, which is built from all relevant fluxes ra" in the system 54 1. 
In the variational approach, the conservation of the individual fluxes; 

Va< = , (1) 

is ensured by means of a pull-back construction based on the notion of a three-dimensional 
matter space. This exercise identifies the spacetime displacements that guarantee ([1]), and 
with respect to which the variation of the Lagrangian is carried out. The detailed procedure 
is discussed in ^0|. For later convenience we simply note that the final result is 

K = n'y.c - - K (v,C + , (2) 

where gab is the spacetime metric and 6gab is the induced variation. 

A key strength of the variational approach is that it correctly identifies the momentum 
/i^ that is conjugate to each flux. This is crucial in a multi-fluid system since the momenta 



should encode the so-called entrainment effect 10|. As an illustration of how this effect arises. 



consider a general isotropic Lagrangian. Taking the view that the fluxes are the fundamental 
variables in the problem, we can build this Lagrangian from the different scalars that we 
can construct. This means that we should consider both 

nl = -« , (3) 

which defines the number density of the x component, and 

nly = -<nl , y ^ X . (4) 

An unconstrained variation of A with respect to the independent vectors and the metric 
Qah then leads to 

M = ^^l5nl + <^'^ ^9ab , (5) 



where the momenta are given by 



with coefficients 



and 



^,l = ga,{B-n'^ + }_^A-^n':^] , (6) 

•4'" = -4"' = -^ . ^^y- (8) 

The momenta are dynamically, and thermodynamically, conjugate to their respective number 
density fluxes, and their magnitudes are the chemical potentials (as we will see later). The 
A^^ coefficients represent the fact that each fluid momentum /i^ may, in general, be given 
by a linear combination of the individual currents ra^. That is, the current and momentum 
for a particular fluid do not have to be parallel. This is the entrainment effect. In the limit 
where all the currents are parallel, e.g. when the fluids are comoving, —A corresponds to 
the local thermodynamic energy density, but in the general case this is not so. 



In terms of the constrained Lagrangian displacements, a variation of A yields 55 1 



(9) 



where we have defined 



fl = 2n^V[,/x^] , (10) 

(and the square brackets indicate anti-symmetrization, as usual). It follows immediately that 
the equations of motion for the individual fluids are expressed as an integrability condition 
on the vorticity (associated with the momentum not the flux!); 

= . (11) 

In we have also introduced the generalized pressure defined by 

X 

Finally, we want to account for the coupling between the matter flow and the dynamics of 



spacetime 



56|. The coupling to gravity follows readily from the fact that the stress-energy 



tensor is obtained as the variation of the matter Lagrangian with respect to the spacetime 



metric. Basically, we know that the geometry side of the problem is obtained from the 
Einstein-Hilbert action, expressed in terms of the Ricci scalar R, 



'EH 



i?v^ d^x . (13) 
Following the standard procedure, this leads to 



5/eh = J GabSg^'V^ d'x = j (^Rab - \gabB^ Sg'^'V^ d'x , (14) 
where Gab is the Einstein tensor. Now, in the coupled matter-gravity system we have 

/ = /eh + /m= [ (^R + A^V^d'^x , (15) 



2k 

where the coupling constant k(= SvrG/c'^) is determined from the correspondence with New- 
tonian gravity in the appropriate limit. This system leads to the usual Einstein field equa- 
tions 

Gab = i^Tab , (16) 

provided that 



2^6 A) 

-g Sg 

or, equivalently. 



^ —I, — , (17) 



T- = ^^iv^. (18) 

Returning to the fluid problem, we see from that the multi-fluid stress-energy tensor 
takes the form 



It is worth noting that when the set of fluid equations, ([T]) and (ITTi) . is satisfied then it is 
automatically true that VaT^ = 0. 

Provided we are given the appropriate matter Lagrangian (a far from trivial problem as a 
realistic model should build on microphysics including the various relevant interactions) we 
now have all the equations we need to describe the dynamics of the fluid system, its effect 
on the gravitational field and vice versa. 

As an aside, it is worth noting that the variational model is more general than the typical 
multi-component models considered in the literature (especially in cosmology) as they tend 



to assume the existence of partial pressures (see 27| for a relevant discussion in the present 
context). 



B. Electromagnetism 



Let us now consider electromagnetism in Einstein's theory. As usual 8| , we construct the 
relativistic version of Maxwell's equations by means of a variational argument with respect 
to the vector potential A"". The corresponding Lagrangian is built from the anti- symmetric 
Faraday tensor; 

F,, = 2V[„A] . (20) 
We also need to couple the electromagnetic field to the matter flow, represented by the charge 



current j". Letting the relevant coupling constant be /^o, the action takes the form |57l | 



'EM 



LemV^ d X , (21) 



with 
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l^EM = --—FabF'''+fAa . (22) 

4/Xo 

However, the current term in this expression is not gauge-invariant. Under a gauge transfor- 
mation of the vector potential, i.e. exercising the freedom to add the gradient of an arbitrary 
scalar field ip, the second term in ( l22l) transforms as 



fA^ ^ fA, + fVa^ = fA, + Va ii^f) - ^ (Vaf) • (23) 

The second term on the right-hand side will contribute a surface term to the action integral, 
and hence can be "ignored" in the usual way. The third term is different. In order to ensure 
that the action is gauge-invariant, we must demand that the current is conserved, i.e. that 

Vaf = . (24) 

The field equations that we derive require that this constraint be satisfied. 

With an action in hand it is straightforward to work out the variation with respect to 
the vector potential (keeping fixed!), and we arrive at the standard result; 

VbF'^" = fiof . (25) 

The relativistic Maxwell equations are completed by 

V[aFbc] = , (26) 

which is automatically satisfied given the anti-symmetry of Fab- 



Finally, a variation with respect to the metric leads to the electromagnetic stress-energy 
tensor being given by 



rpEM _ 1 



g^'^FacFbd — -gab {FcdF'^'^) 



(27) 



It is worth noting that this leads to 

V^T^t = JaF"' = -ft , (28) 

which (as we will see later) defines the Lorentz force /£. 

In principle, the electromagnetic dynamics is now fully specified, as we can solve the 
system for the vector potential A"-. However, in most applications it is more intuitive to 
work with the electric and magnetic fields E°- and 5". The down-side to this is that these are 
observer dependent quantities. This is obvious since varying electric fields generate magnetic 
fields and vice versa, and the induced variation depends on the motion of the observer. 

According to an observer moving with four-velocity m", the Faraday tensor can be ex- 



pressed as [58 1 



Fab = 2U(^aEb) + eabcdU'B^ , (29) 

(where round brackets indicate symmetrization). This defines the electric and magnetic 
fields as 

Fa = -u'Fba , (30) 

and 

The physical fields are both orthogonal to u"", so each field has three components, just as 
in non-relativistic physics. We also need an expression for the current, and it is natural to 
decompose this in a similar way; 

f = au" + r , where J"u„ = . (32) 



C. A comfortable marriage 

So far, we have done quite a lot of preparatory work, going over standard territory without 
adding any real new insight. Our patience with this exercise is about to pay off, as we will 
now be able to make swift progress. This illustrates the advantage of having a well-grounded 
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action principle for coupled fluids, and an identification of the true momenta, and shows 



how easy it is to incorporate electromagnetism into the multi- fluid system [28|]. We simply 
need to consider multiple charge carriers with identifiable fluxes, n^, and individual charges, 
gx, such that the charge current associated with each flow is 

Jx = ?x< , (33) 
and the total current, that sources the electromagnetic field, is given by the sum 

X 

It is worth recalling that the variational derivation in Section JIB requires that the current 
is conserved. However, this constraint is automatically satisfied if each individual current 
is conserved, as assumed in the variational multi-fluid model. Hence, we simply have to 
change the electromagnetic Lagrangian to 



L 



EM 



+ (35) 



to combine the two models. 

It is easy to see that the equations that govern the electromagnetic field remain exactly 
as before. However, the coupling to the current leads to modified fluid momenta; 

^l = ^^l + q'^Aa, (36) 

which satisfy the equations of motion 

2<V[,/i^] = . (37) 

As an alternative, we can write this as an explicit force-balance relation. Moving the elec- 
tromagnetic contribution to the right-hand side, we get 

f: = 2n^V[,/i^] = q^nlF,, = jIF,, . (38) 

To see that this result makes sense, note that the total energy-momentum tensor is easily 
obtained as the sum of the two previous expressions; 

= < + T^l . (39) 
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This means that we must have 



ab 



V T 

Va-i M 



Va-^EM 



In other words, the combined system is such that 



fb — ' 



(40) 



(41) 



The variational formahsm naturally lends itself to a consideration of conserved quantities, 
like the magnetic helicity j9|. The discussion becomes particularly elegant if carried out using 



the language of differential forms 



591]. We will not discuss conservation laws in this paper. 
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30| 



but the interested reader can find relevant recent discussions in 

Before we proceed, it is worth digressing on the fact that the charge current does not 
enter the electromagnetic stress-energy tensor fl27|) . As this is a key (albeit somewhat tech- 
nical) point, it is worth demonstrating the result in detail. To do this, let us focus on the 
contribution to the total action from the matter-field coupling; 



Variation of the integrand then leads to a sum of terms of the form 



(42) 



(43) 



Naively, the first term affects the Euler equation, the second leads to the current term in 
the Maxwell equations and the final term should enter the stress energy tensor. However, 
the last contribution is cancelled by a term originating from the variation of the matter flux. 
Using ([2]) in fH5]) (ignoring surface terms) we arrive at 

6 «A,v^) = (2e>''V[,A] + K6A,) . 



(44) 

The first term enters the Euler equations ( l38ll and the second leads to the current term 
in the Maxwell equations. The electromagnetic contribution to the stress-energy tensor is 
completely determined by the first term in fl22p . leading to f l27p . 



III. RESISTIVE MAGNETOHYDRODYNAMICS 



The formalism developed in the previous section provides a general framework for de- 
scribing the dynamics of charged multi-fluid systems in relativity. However, as the model 
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arises from a variational analysis it does not account for dissipative mechanisms. Hence, we 
need to amend it if we want to model, for example, resistivity. This is obviously a key aspect 
if we want to be able to model the evolution of electromagnetic fields in various astrophys- 
ical and cosmological settings. However, we know that the general dissipative problem is a 
severe challenge in relativity. We also know that many different dissipation channels may 



affect a generic multi-fluid system [23|, l32|. Hence, we set a more modest target and explore 



the role of a simple, phenomeno logical, resistivity. As it turns out, the problem involves 
tricky issues already at this level. In general, any dissipative mechanism will generate heat. 



so a realistic model must account for the associated heat flux. However, t 
known to be associated with both causality and stability issues in relativity 



lis problem is 
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problems can be resolved 



21 



19|. These 



22|, but we must proceed carefully. 



Given the various issues involved, we consider the resistive problem at two different 
levels of sophistication. First (in this section) we consider a cold plasma, where the various 
relative velocities are sufficiently low that the problem simplifies. Having understood this 
problem we proceed (in the next section) to consider the general problem, with arbitrary 
velocities and the presence of a heat flow. In each case, we consider a system with two 
charge carriers, with individual particles carrying a single unit of charge. This means that 
the models apply to both pair-plasmas with positrons and electrons and proton-electron 
plasmas. These examples provide useful illustrations and the discussion serves to highlight 
differences between the two systems. 

Throughout the discussion, we assume that the system is fully ionized. That is, we do not 
allow for the presence of a charge neutral component, as would be required if we wanted to 
model a magnetized neutron star core, for example. The inclusion of such a component is, 
in principle, straightforward although the algebra obviously gets more involved (especially 
if one accounts for entrainment). 



A. Choice of frame 



We focus on a two-component system with one component, labelled p, carrying a single 
positive unit of charge q'p = e while the other component, labelled e, carries a single negative 
unit of charge qe = —e. The associated charge currents are = en'^ and = — era^, 
respectively. We will not, initially, make any assumptions regarding the relative masses of 
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the two components. This means that the model apphes to both pair plasmas and proton- 
electron systems (indeed, any two-component system with electrons and single charged ions). 

A key aim of the exercise we are embarking on is to derive the relativistic version of Ohm's 
law. Basically, we want to start from the charged two-fluid system and arrive at a model 
from which the assumptions associated with standard relativistic magnetohydrodynamics 
become clear. This discussion will obviously involve both the electric and the magnetic 
field, as well as the charge currents. Now, we know that E"" and B"^ are observer dependent 
quantities. Hence, the model involves a judicious choice of observer. It is natural to begin 
by considering this issue. 

Given an observer with four velocity u"" (normalised such that u'^Ua = —1) we can de- 
compose the two fluxes; 

= tIxm" , where nl = -n'^nl and u^ul = -1 , (45) 

using 

< = 7x (m" + <) , where mX = , and 7x = (l - vl)'^^^ . (46) 

In the first instance, we will assume that the "drift" velocities v'^ are small enough that we 
can linearise the model, i.e. assume that 7x ~ 1. This model should be relevant for cold 



plasmas [60 1 



We also need the fluid momenta, which would generally involve entrainment between the 
two components. However, as we are not aware of a physical argument for the presence 
of entrainment between protons and electrons (or, indeed, positrons and electrons), we do 
not account for this effect here (although we will consider it when we discuss heat flux and 
entropy later). This means that we have 

= B^nl = B^^.n^ {ua + <) . (47) 

The chemical potential of each component is generally defined by 

/xx = = nxi3^ , (48) 

which means that, 

^4 = /ix7x {Ua + O ^ ^x {Ua + <) • (49) 
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With these definitions we can write the (hnearised) fluid stress-energy tensor as 



= -^Qab + (np/ip + rie/ic) UaUb + (rip/ipi;^ + U^fleVl) Ua + {u^^-pVl + U^fleVl) Ub ■ (50) 

Contracting this with u"" we get an expression for the momentum flux; 

U^T^ = - Upfip - nefJ,e) Ub - {npflpV^ + riefievl) . (51) 

Contracting with again, we flnd that the energy density measured by the observer is 

p = u^u'T^ = -^ + ripfip + ne/ie . (52) 

We see that, in the hnear model \1/ is the pressure. Hence, we replace it with P in the 
following, leading to the anticipated thermodynamic relation (the integrated flrst law) 

P + P = ripfip + UePe ■ (53) 

Note that, as we are only considering the fluid contribution here, our deflnitions of P and p 
do not include electromagnetic effects (i.e. the magnetic pressure is not accounted for yet). 

From flSTl) we see that we can choose observers such that there is no relative (fluid) 
momentum flux by setting |6l| 

npPpVf + UePeVl = . (54) 

This leads us to deflne a velocity such that 

(P + p)v'' = UpPpV^ + riePeVe , (55) 

and highlights the relevance of the frame in which = 0. We express the second degree of 
freedom in terms of the relative velocity 



With these deflnitions we have 



and 



= < - . (56) 



f - , (58) 
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and the charge current takes the form 

T) T) 

r = e (np - ne) (n'^ + t;") + e-^^ (/ip + p,) w'' . (59) 

From this result we read off the charge density a = e (rzp — Ue) in the observer's frame. If 
we assume that the system is charge neutral on macroscopic scales, a natural assertion for 
systems where the charge carriers (like the electrons) are highly mobile and one of the key 
assumption in standard magnetohydrodynamics, then the current simplifies to 

f = r = e^^{p^ + p,)w^. (60) 

Moreover, in the case of a charge neutral plasma we have P + p = (/ip + Pe) which 
means that the current takes the final form 

r = euew" . (61) 



B. The resistivity 

In order to account for the resistivity, we need to add a phenomenological "force" term 
to (1381) . This additional term should represent the dissipative interaction between the two 



components, and from non-relativistic intuition [11, 12|, we expect it to be linear in the 
relative velocity between the two components. We also see from (138|) that the required force 
must be orthogonal to each respective flux (note that this condition must be relaxed if we 
want to allow for particle creation/destruction). Based on these points, we let the resistive 
forces take the form 

f^ = en±fnl = -n±'^'jt, (62) 

and 

/: = en nl = n ±f J, . (63) 

These expressions represent linear scattering of the two components. The resistivity expe- 
rienced by one component is proportional to the number of particles of the other kind that 
flows relative to it. 

The resistivity is further constrained by the fact that the sum of the forces must vanish 
(essentially Newton's third law). This follows immediately from the fact that the divergence 
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of the non-dissipative stress-energy tensor [which arises from the sum of (l38l) ] must vanish. 
With the suggested forces, we have 

/; + /: = en {l.f nl+ nX) ^ eU {n^ - n,) . (64) 

This shows that the hnearized model is only consistent as long as the system is charge 
neutral. If there is charge imbalance, we need to alter the model. At first sight, this 
may seem surprising but it is actually quite natural. The model only accounts for the two 
charged components, whereas the general system would also have the heat generated by the 
dissipation. The correct interpretation of (IMI) is that, for a charge neutral system, there is 
no heat generated at the linear level. In order to consider a more general system, we need 
to account for the heat. Then the force balance is ensured by introducing an additional 
component, which we will take to be the entropy, with a corresponding force of the required 
form. We will discuss this extended system in the next section. For now, we simply assume 
charge neutrahty and note that the corresponding low-velocity model describes a "cold 
plasma" in the sense that there is no heat generated in the system. 



C. Generalized Ohm's law 

The problem under consideration has two fluid degrees of freedom, represented by (!38|) 
with the added resistivity terms (on the right-hand side). One can (obviously) combine 
these two equations in different ways. It seems natural to adapt the standard strategy from 



non-relativistic plasma physics (ill . Il2l ] and consider a "total momentum" equation alongside 
a suitably weighted difference. The first of these equations follows by adding (155]) . and from 
the discussion in the previous section we know that this leads to 



VaT''' = , (65) 

as the sum of the resistive forces vanishes (to linear order). In order to represent the 
second degree of freedom, we divide the two equations from fl38|l by n^fix and then take 
the difference. The weighting (different from that used in other recent discussions of the 



problem [l5j) is motivated by the Newtonian limit, where fix fn^ (the rest mass) and 
corresponds to the " centre-of-mass" frame. With this weighting the difference equation 
simplifies considerably. One may obtain the same final result with a different weighting, but 
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the analysis would then have to make explicit use of the total momentum equation (l65l) in 
simplifying the expressions. Our route is more direct. 

The difference equation that we require is made up of three pieces. Considering first the 
fluid contribution to fl5S|) and the definition of the chemical potentials, we have 

2<V[a/U6] = nx Va/ix + nx/ix<Va< . (66) 

We also have 

±f = g^^ + ulul +2u^%^^ , (67) 

where l_"-^= g"-^ + u°'vf' is the projection orthogonal to u"". Using these results, we find that 
the weighted difference (let us call it /^) in the linear model takes the form; 

^ U^VaWb + w'^VaUb + — ±pfe Va/ip - — Kb Va/^e • (68) 

The last two terms expand to 



" +u''Vb + Ufot;'') ( — Va/ip - — VayUc 
+ _^ ^ (^"^fc + "fc^") K/^gVa/ip + np/ipVa/ie] • (69) 

This expression obviously simplifies somewhat in the frame where v'^ = 0. The result 
also simplifies for the two specific examples we are considering. For pair plasmas we have 
yUp = /io = yU so in the chosen frame the final result would be 

2 

fF = WVaWb + W^VaUb + -n(aU7b)V> . (70) 

Meanwhile, for a proton-electron plasma we may assume that /Xg ^ /^p in which case (when 



= 0) we are left with 62 1 

fF ^ U^'VaWb + w"VaM6+ l-\ (— Va/ip - — Va/ig ) + — M(aW;b)V>e • (71) 



In each case we can replace the relative velocity w"" with the charge current via Eq. ([6T 
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The expressions we have obtained represent the left-hand side of the equation that we 
are developing. The right-hand side is made up of two pieces. The first is the weighted 
difference of the two magnetic forces from fl38p . That is, we have 

= e(^-^)(E. + ,fF^)-t^^4^yfF^. (72) 

As before, this simplifies in the frame where v"" = 0. In addition, for (charge neutral) pair 
plasmas the second term vanishes identically and we are left with 

/r = e^^^^E, = ^E, . (73) 

/ip/ie /i 

Meanwhile, for a proton-electron system we would have 

^-{E,- w^F,a) . (74) 

In this case we need to consider the remaining term involving the Faraday tensor in more 
detail. From (129|) it is easy to see that we will have 

w'^Fba = Ub (w'^Ea) + ebcdw'B'' , (75) 

where we have defined ebcd = ^nhr du"' in order to make the final term resemble the standard 



three-dimensional cross product 63|. Since w"" is proportional to the charge current, we 



recognize the two terms as the Joule heating and the Hal 
are notably absent in the linear model for a pair plasma 



effect, respectively. These effects 
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Finally, we need the weighted difference between the two resistivities. That is; 

f R _ fP f e _ ^'^ |P a ^'^ le a 

Jb — Jb Jb — -"-afe "'e -"-ab "'p 

^p/^p ^e/^e ^pf^p ^e/^e 

1 ,P , 1 ,c^,a^ ^ P + P 



Vnp/ip ne/ie / (np/ip)(ne/ie) 

where the fact that the current is proportional to w"' in the charge-neutral case has allowed 
us to neglect the various terms in the projections. 
The final result now follows from the combination 

/° = IT + If ■ (77) 
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In the general charge neutral case we have, in the frame where v°' = 0, 



e J^b \Ub[W r^a) + ^bcdW D — K- ^ r-Jfe 

= U^VaWb + w''VaUb+ l-\ ( — V^/ip - — Va/ic 

2 

+ .p X Uf^aWb) [ne/ig V/ip + rip/ipV/ie] • (78) 

This result can be simplified by projecting out the contribution along (which will not 
affect the expression for the electric field). This leads to the "final" result; 

(P + p) e(/ip-/ie) , d P + P 
e J^b €bcdW n — K- — -Jh 



PpPe (rip/ip) (ne/ic 

U^'VcWa + — Va/ip - — Va/ie 
f^p Pe 

rip 

+ 



+ W^'VaUb 
WbU" [/ieVa/ip + /iJVa/ie] ■ (79) 



(P + p)fipfle 

In these expressions, it may be useful to decompose VaWft in the standard way, see for 
example That is, we use 



VaMfo = (Tab + ^ah " UaUb + -0 ^ab , 

in terms of the expansion scalar 
the shear 

(Tab = D{aUb) , 

where the angle brackets indicate symmetrization and trace removal, and 



DaUb 



-a b 



^cUd ■ 



(80) 

(81) 
(82) 

(83) 



64| 



The merit of using this (totally projected) derivative is that the individual terms in f ISOl) are 
perpendicular to u". We have also defined the vorticity 



l^ab = D[aUb] ■ 



(84) 



The decomposition flHOj) makes the coupling between the charge current J"' and the nature 
of the fiuid motion more explicit. 
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The final relation for pair plasmas can now be written 65 1 



F ^ T 



ab J" + J" ( (Tab + (^ab + 77^ -"-ab 



4 
-t 

3 



(85) 



where the dot represents the comoving time derivative u'^Vc- As already mentioned, this 
expression is notable for the absence of the Hall effect, i.e. there is no term proportional to 

^abcJ''B'^- 

The case of a proton-electron plasma is only slightly more complicated. After neglecting 
/ie compared to /ip, we end up with 



Eb e^cd.J'^B'^ 

eup eup 



Jb 



e^Hp 



-ab J" + J" \ (Tab + ^ab + ^ab 



V,/ie. (86) 



In this case, the Hall term is obviously present. We also have a "Biermann battery" term, 
1\ Va/ie, which would serve to generate a magnetic field even if there was no field initially 
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It is easy to show that our final results agree perfectly with the results obtained in 15|. 

Before moving on, it is useful to consider the relation between our results and the common 
starting point for discussions of resistive effects in numerical simulations 3J-|36| . Much of the 
relevant literature builds on the work by Bekenstein and Oron 33| . Ignoring the right-hand 
side of Ohm's law in the proton-electron case we have 



1 IZ IZ f \ 

Eb = ebcdJ'^B H Jb = -Lab +7F;^bacdU^B I J" = SbaJ"' ■ 

erie erip UgC \ l-i 

Define ( = 1/TZ and a = TZ/rice to get 



(87) 



Sba = — (-Lfea +C^bacdU''B'^) . 



Inverting this, we arrive at 



Jb = o'"'^Eb , 



with 



ab 



a 



a 



^^ab ^^2^a^b _ ^g^fecd^^^^) . 



(88) 



(89) 



(90) 



1 + CB^ 

This is the result stated in |33], once we account for the different sign conventions. At this 
point, we can make an important observation. It is easy to identify the Hall effect in the 
initial expression ( 187|) . but its presence is more convoluted in the alternative expression (l90l) . 



That this can lead to conceptual confusion is evidenced by 



37l |. where numerical evolutions 
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for a truncated form of (pOj) are carried out. The considered model includes a peculiarly 
amputated Hall effect, the actual meaning of which is unclear. This lesson tells us that an 
understanding of the physical origin of the model is imperative. 



D. Towards ideal MHD 

We are now in a position where we can assess the relative importance of the different terms 
in the generalized version of Ohm's law fl79p . Let us first consider under what conditions we 
can neglect the inertia of the charge current compared to the resistivity. In order to do this 
we need 

J>-^J^^J. (91 

It is natural [33'] to associate the resistivity with a relaxation timescale such that 

7^ = ^ . (92) 

If we also assume that the the dynamics has a characteristic timescale r^, such that J ~ ^/t^, 
then it is easy to see that we can neglect the inertia of the charge current as long as 

Td > . (93) 

When this condition holds, i.e. for sufficiently slow dynamics, the system is essentially obliv- 
ious of its plasma physics origins. This condition shows why ideal magnetohydrodynamics 
is a good model for slowly evolving, or stationary, systems. 

Next let us compare the Hall term to the resistivity. The former dominates (in magnitude) 

if 

5 > 7^ , (94) 
which, if we introduce the electron cyclotron frequency 

O0c = — , 95 

leads to the condition 

UcTr > 1 . (96) 

This means that the electron executes many cyclotron "oscillations" before the motion is 
damped. 
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Finally, we need to establish when the resistivity can be neglected compared to the electric 
field. This requires 

E^—J = -^J. (97) 

Here we need to make use of Maxwell's equations (see below), which lead to (in Gaussian 
units!) 

E ~ ^fi , (98) 

where we assume that the dynamics has a characteristic length-scale L, and an associated 
velocity V = L/ra. This leads to the final condition 

J <ti rieC ^ — ^ (jJcTr <^ ■ (^9) 

The last condition is required since we also want to be able to neglect the Hall term. We 
are essentially left with a low- velocity constraint. 

These rough estimates provide useful insight into the applicability of "ideal" magnetohy- 
drodynamics, which corresponds to the assumption that E"" Q. The usual argument for 
this is that the medium is a perfect conductor, i.e. 7^ — 0. However, this limit only affects 
the resistive term in (1791) . We still have to argue that the remaining terms are unimportant. 
This is not quite as easy. At the end of the day, ideal magnetohydrodynamics is more an 
assumption than an approximation [ll] (the interested reader may want to compare the 



present discussion to the variational derivation of magnetohydrodynamics in [38l|). 



E. The remaining fluid equation 

So far, we have focused on the weighted difference between the two momentum equations 
in the plasma. To complete the "single-fluid" model we need to express the remaining degree 
of freedom in terms of our chosen variables. It is natural to obtain the required equation 
from (1651) . 

As a flrst step, we consider the non-magnetic contributions. As the individual number 
fluxes are conserved in the variational approach, we see that 

VaP +iP + p)u' + 2u^''v'JnpVafip + 2u^''vSvy alie 

+ n^^i^{vl + vlVay')+n,^i,{vl + vtVav') . (100) 
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The first line is exact, while the second line holds at the level of linearised relative velocities. 
Expressing this result in terms of and w"", we have 

+ '^^u'w'^(-VaP^--VaP^+{P + p){v' + v'^Vau') . (101) 
1 -r p y /ip Pq J 

Here we have used the fact that we are considering a charge neutral system. This result 
provides the left-hand side of the final equation. As discussed in section IIC the matter 
contribution is balanced by the electromagnetic stresses, which provides the right-hand side 
for the equation we are interested in. This takes the form 

- V„T^^ = a^'^) + e^'^'jaB, . (102) 

The final equation will only have spatial components with respect to m", but the relations 
we have written down so far also have parallel components. However, if we contract the 
combined equation with Uh and compare to what we get if we contract our generalized 
Ohm's law (1791) with the current, then we see that the two results agree (at the linearized 
level, of course). Hence, only the orthogonal component contains new information. In the 
frame where f " = the final fiuid equation takes the form 

(P + p) L'^'' VaP = e'^'JaB, . (103) 

This is simply the perfect fluid equation of motion augmented by the Lorentz force. 

To complete the model, we may also consider the two conservation laws ([1]). It is straight- 
forward to show that the difference between these corresponds to the required conservation 
law for the charge current. Meanwhile, after making use of the component aligned with 
from the weighted difference equation that leads to ( 1791) , the sum of the two conservation 
laws can be written 

J^E^ = ^±P{n^ + n,e) . (104) 

He 

Moreover, one can show that (at the linear level) this expression also follows from the 
component of the total momentum equation fl65|) that is aligned with m". 

At the end of the day we have two scalar equations and two equations governing velocity 
components that are spatial with respect to u"". Thus, we have explicitly accounted for the 
degrees of freedom of the original two-fluid system. To complete the system we also need 
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Maxwell's equations. Before discussing these, let us make a brief diversion and touch upon 
a model that is common for neutron star magnetospheres. 

The conditions in the magnetosphere of a neutron star (or, indeed, a black hole) are 
expected to be such that there is sufficient plasma present to carry a charge current, but 



the associated inertia can be neglected [39|, |40|. Thus, we may neglect the inertia in f llU3l) . 



essentially decoupling the matter from the magnetic problem. This leads to what is known 



as force-free electrodynamics 
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42|. In these circumstances we would have 

e'^V.^e ^ , (105) 



which implies that charges may only flow along the magnetic fleld. 

The force-free assumption can obviously be used independently of the assumptions that 
lead to ideal magnetohydrodynamics. Basically, one may envisage a range of different "ap- 
proximations" depending on the circumstances. The force-free model simplifles magneto- 
sphere modelling, but one must apply it with care since it breaks down near magnetic 
neutral points. In the context of the present discussion, it is also worth noting that fllOSp 
may be extended to include various dissipation channels (like shear viscosity) other than the 
pure coUisional resistivity that we have accounted for. If the multi-fluid aspects are taken 



seriously 25|, this may lead to a much more complex problem. 



F. Maxwell's equations 

Given an observer moving with m", representing a flbration of spacetime, the decomposi- 
tion of Maxwell's equations is standard. Nevertheless, we list the results here for complete- 
ness. For a more detailed discussion, see for example, jsl. 

First of all, 

V„F^'^ = /ioi' , (106) 

leads to 

WbEa = VaE'' - E'^ila = flQCT + e'^'^UabBc = f^oCT + 2W''Ba , (107) 

where we have deflned the vorticity vector as 

W'' = ^e'^'^cohc , so that cOab = eabcW , and u'^Wa = . (108) 
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We also get 



Secondly, 



leads to 



and 



-ab E' - eabcV'B' + ^loJa = ( - Uab - ±ab ) E' + eabcU' . (109) 



V[,F,,] = , (110) 



V,B, = -2W''Ea , (111) 



^ab B' + eabc^'E' = -eabcU'E' + (^aab - COab - ^0 ±abj B' . (112) 

It is easy to see that, if we consider an inertial observer, these results reduce to the 
standard text-book form of Maxwell's equations. The complete expressions given here are, 
of course, useful if we are interested in more general settings. In particular, they highlight 
the coupling between the electromagnetic field and a given fluid flow (with shear, vorticity 
and expansion). 

IV. ADDING ENTROPY: HOT PLASMAS 

The low-velocity model we have discussed so far is consistent and applicable to many 
situations of interest. It also provides a number of potentially important extensions of the 
ideal magnetohydrodynamics that tends to be used in relativistic astrophysics. However, as 
we have already hinted at, the model does not account for the presence of heat. This is an 
unfortunate omission since resistivity is a dissipative process and hence will be associated 
with entropy variations constrained by the second law of thermodynamics. This effect turns 
out to be quadratic in the relative velocities, which is why we got away with neglecting it in 
the linear model. In a more general setting, we need to account for the induced heat flow. 
The problem of heat in relativity is, however, known to be thorny. A model needs to be 



constructed carefully in order to avoid unwanted instabilities and causality violation [16Hl9|. 



As recently demonstrated, one can construct a satisfactory model by treating the entropy as 



an additional fluid component 66| , and accounting for entrainment between the entropy and 



the other components in the system 2l|, |22l]. This entropy entrainment is closely associated 



with the inertia of heat and the flnite thermal relaxation timescale that is required in order 
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to avoid superluminal signal propagation. We develop our model with these key points in 
mind. 



A. Setting the stage: A three-fluid system 

We consider a hot plasma consisting of the two charged components from the cold model, 
labelled p and e as before, and an additional entropy, which we label s. As the entropy plays 
a special role, being constrained by the second law, we will single out this component by 
letting its flux be given by = while the corresponding chemical potential is Ga = /i^. 
The latter determines the temperature measured by a given observer. With these definitions 
we have the total stress-energy tensor {l^ 

Tah = ^9ab + + nlnl + SaQb , (113) 

where the generalized pressure \& is defined as 

^ = A-</iP-</^:-s'^0a • (114) 



Following [2l|, |22|] we account for entrainment between the entropy and each of the material 
components (encoded in coefficients A^"^). Thus, we have the momenta 

/i^ = SX + ^^^Sa, x = p,e, (115) 

and 

e, = i3^s„ + ^pX + ^'X • (116) 

As in the low-velocity model, we introduce a family of observers that allow us to define 
the electric and magnetic field components. We now have a number of different options. 
The strategy that we adopt provides a natural extension of the Eckart frame for a single 



component matter model, cf., [21|, 



. To be specific, we choose the observer frame to be 
such that the only relative momentum flow is due to the heat. 

Defining first of all the number densities as measured in the respective fluid frames, we 
have [63] 

nl = -KK, and = -s'^s^ . (117) 

Decomposing the velocities with respect to a specific observer moving with -u" we then have 

< = n.7x {u' + <) , uV, = , ^,= {l-vl) , (118) 
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leading to the number density measured by the observer being given by 

nx = -Ma< = nx7x • (119) 

Similarly, we have 

= (w'^ + <) , u%: = , 7s = (1 - v!) '^'^ , (120) 

and 

S = -Uatt" = s% . (121) 
It is also natural to introduce the chemical potentials inferred by the observer; 

/ix = -M>a = ^xi?" + SA""' , (122) 

and 

Q = sl3' + ripAP' + n.A''' , (123) 

With these definitions it is straightforward to show that the total energy density measured 
by the observer will be 

p = u^u^Tab = + rip/ip + ne/ie + sQ , (124) 

corresponding to the (integrated) first law of thermodynamics once we identify \E' as the 
generalized pressure and as the temperature. Meanwhile, the momentum fiux relative to 
the observer's frame is given by 

u^Tah = -pub - ripPpvf - UePeVl - sQvl . (125) 



Here it is, first of all, natural to identify the heat fiux as 2l|, |22| 



g'^ = sOv^ . (126) 

We also see that we can choose the observer frame in such a way that this is the only relative 
momentum fiux. To do this, we let 

(p + ^) = n^ppv; + ne/ie< = . (127) 

This is the natural extension of the "centre of mass" frame we used in the low-velocity 
model, cf., eq. fl55|) . We also define the velocity difference (as before); 



(128) 
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In the frame where f ° = (which will be assumed from now on) we have 



< = , (129) 



and 



""P^P w'^ , (130) 



which means that the charge current can be written 

r = e (np - ne) + e-^^ (/ip + ft,) w'' . (131) 

At this point we recognize an important difference with respect to the low- velocity discussion. 
While we were naturally led to the assumption of charge neutrality in that case, the situation 
is much less clear now. This is immediately obvious from fll3ip once we recall that the 
densities in the first term are measured by the chosen observer, not in the respective rest 
frames. Hence, it would not be appropriate to assume that rip = n^. at this point. 

B. Friction and causal heat flow 

Having introduced the various ingredients, let us move on to the new aspect of the 
problem; the equation that governs the heat propagation. As we are treating the entropy as 



nerma. 
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dynamics will 



22j that this will 



an additional fluid, it follows from the general analysis in [l^ that the t 
be governed by its own momentum equation. We already know from 
lead to an equation that contains the thermal relaxation that is required to ensure causality. 
However, as the entropy need not be conserved this momentum equation takes a slightly 
different form from those that govern the (individually conserved) material components. We 
have |2l|,|22| 

2s''V[„efe] + e^rs = n , (132) 

where 

Vas" = r, > , (133) 

in accordance with the second law. 

Building on the analysis of the low-velocity case, we know that the overall conservation 
of energy and momentum requires 

E/x = o f: = -f^-f:- (134) 
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The form of the "force" that acts on the entropy thus follows immediately from the forces on 
the charged components. Extending the low-velocity model, we will allow for resistivity due 
to scattering between both charged species (7^) and entropy (5x). Thus, we let the forces 
take the (still phenomenological) form 

f; =l.f {eUnl + S^ul) , (135) 

and 

/: =±f (e7^< + S^ul) , (136) 

where we have used s'^ = sm". Combining these expressions and expanding the projections, 
we arrive at 

f: = -en [n, - np7e (l - + <) - eU [n^ - n^^l (l - v\vt)\ (w" + t;^ 

- 5p {ul - 7.7p«; (1 - v\vt)\ - [< - 7s7e< (1 - v\vt)\ • (137) 

Let us now return to fll32p . focussing on the entropy creation rate. Contracting the 
equation with the observer's four velocity we easily arrive at 

er^ = -u^'fi + isu^y^^^Qv, • (iss) 

We need to constrain this in such a way that the right-hand side is non-negative. To do this, 
we first need the contraction between the entropy force and the four velocity. This leads to 

6 (2 2,2 2\ 2 I 7s Tc 2 / \2 , c 2 / n21 2 

Jb = ("e/^e7p + ^p/^p7e j ^ + ^ ^^2 ['^P% (^e/^e) + S^-f^ (Wp/ip) J W 

+ -^^^ ('5e7e"'p/^P - 5p7jne/ie) • (139) 

In this expression, the first two terms on the right-hand side will be positive as long as 7^ > 
and iSx > 0. The sign of the third term is not so clear. 

Moving on to the final term in (11381) . we first of all note that it is proportional to cx q"'. 
Defining 



and 



77 77 

^2 = ^ {^^eA^' - f^,An , (141) 
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we have 



and we find that 



(142) 



2«''V[,e,] = - ±^ - - - (32Wa - PlQa - $2Wa ' (Ag' + V . (143) 

Combining ( I139p and ( 11431) we see that Fg satisfies the required constraint provided 



that 
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-K 



7s 



s(p+ ^) 



(5e7fnp/ip - 5p7pne/ie) U)a , (144) 



with K > 0. This has the form of a Cattaneo-type equation 2l|, and a comparison of the 
and Qa term suggests that the thermal relaxation time is 

kBi 

T = -V . 145 

However, the equation is also coupled to the four-acceleration -u" and the variation of the 
charge current, in terms of w"", so if we want to infer the actual relaxation times in the 
problem we need to consider the coupled system. 

Combining the relevant contributions, we find that the total entropy creation rate is given 

by 



^ + (^eVe7p + np/ip7e) 



en 



+ 



Is 



> . (146) 



Ohm's law 



The derivation of the generalized form of Ohm's law follows the same steps as in the 
linear model, although now we need to keep careful track of the different redshift factors, 
etcetera. Basically, we want to construct the weighted difference between the momentum 
equations for the two charged components, but the two momenta now depend also on the 
entropy flux. In the frame associated with our chosen observer, we have 



(147) 
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where we have introduced the convenient combinations 



sQ p + ^ 



Wr 



and 



sB p + w 



(148) 



(149) 



(obviously expressed in the chosen frame). 
Given these expressions it follows that 

+ ±t Va/ix + n^W'vfVali. + 2n,M'^V[aW,^ + 2nXV[a/ib] • (150) 

The weighted difference equation (inevitably rather complicated) combines three pieces. 
On the left-hand side we have 

+ ±t ( — Va/ip - — Va/ie ) + ( — t;pVa/ip - —V^V af-ie 
V/ip /ie / V/^P f^e 

+ -w'^Vi.WP - -«"V[„>V,^] + -t;p"V[,p5 - -<V[,/i^] . (151) 

/^p /^e /^p /^e 

In the frame associated with the chosen observer, this expression takes the form 

/fe^ = Wb + w"VaMb+ ±b ( — VaPp - — VaPe 

H ^ — — rv2M(aWi,) [ne/ig V/ip - rip/ipV/ie] 



69| 



/ipPe(P + ^)' 



P+^J Pp 



-V>p - 



p+^y Pe 



2Pt."V[„^5] - 2^«;'^«;[,V.] ( ^ 
p + W \p + w 



+ 2 



ricPe 



p + ^ 



^ppp 
p+ ^ 



O ^ppp „,,av7 -lAio 



(152) 



Here we have introduced 
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Pe/ip(p+ ^) 



(neP,MP^ + nppJ^-) 



npTie 



e 



p + ^ 



?T-ePe ^pPp n 

-LJ H P2 

^ppp^cpe 



(153) 



30 



and 



V 



We have also used the projection orthogonal to w"; 



(154) 



(155) 



±;ffe= W gab - WaWb . 



(156) 



Meanwhile, the right-hand side is made up of, first of all, the combined friction forces; 



1 



np/ip 

en 



UeUe 



e 
-ah 



p + 

7s7r 



/S22 '\ 2 1\ 2 /2 2 2 2N i -i 

^ (^p/^p7e - ^e/^e7pj ^ Mfo - (n^/ie7p + ^p/^p7c j ( 1 



w gg _ ricHcW^ 
sQ p + ^ 



Ub + 



1 - 



W 



P + ^ 



1 ^ — I 



p + ^ sG 

2~ 



nePe P+^'V'50 P + ^ 



1 - 



_rapPp_ 
p+ ^ 



to 



sQ 



^p/ip 
p+ ^ 



nppp 



where the second equality holds in the frame where v"" = 0. 

The final part accounts for the electromagnetic field. We need 

ripPp ^ riePe yWp/^e V/^P /^e / 



ba 



Wb 



(157) 



(158) 



/ip/ie / /ip/ic (P + ^) 

Again the second equality holds only in the chosen frame. 

The final relation follows from the combination ( 1771) after projecting out the component 
orthogonal to u"". The result is (inevitably) rather complex and may not be particularly 
instructive. Yet, we provide it in the interest of completeness. The generalised form for 
Ohm's law for a hot two-component plasma can be written (expressed in terms of w"" rather 
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than the spatial component of the charge current J", for convenience) 

Wb + w"Va'U6+ ±t ( — Va/Xp - — Va/^e ) 



— Va/ie 



e7^ 



(^e/^e7p + ^p/^p7e) ( 1 ^p^^" ) ^6 



np/Xp 



1 - 



1 - 



n 



where we have defined 



V /^p/^e 



Eb + e 



-^6,g,^«S% (159) 

(160) 



D. The total momentum equation 



As in the two-component system, the model is completed by the total momentum equa- 
tion, which follows (more or less) immediately from the divergence of the stress-energy 
tensor. In the frame moving with u"' (where — 0), the matter stress-energy tensor takes 
the form 

= pUaUb+ ±ab * + 2'U(ag6) + aWaWb + '^^(aQb) + ^QaQb , 

with 



fig/igfipPp . 

a — ^z-^ 1 



HpPp 



;i6i) 



^ n^PeripPp A _ ^pie\ _ s (ne//e - ripPp) 
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As usual, VaT"'' = can be divided into a component along the four velocity and an 
orthogonal piece. After a bit of algebra, recalling that the electromagnetic contribution is 
given by eq. (125]) . we find that the former can be written 



p + (p + ^) VaW" + Vag" 
- Ub 



q' + aw^Vaw' + ^ [w^Vaq" + q^V aw') + ^q^V aq' 







ripne 



p + ^ 

Meanwhile, the orthogonal projection leads to the momentum equation 

(p + ^) L""^ Va^ + g^VaM^ + q^Va + V^) 



(Pp + Pe)K^") . (163) 



e (rip - rie) i?" + (pp + pe) e^'^^u;^!?, , (164) 



p + ^ 



where we have defined 



and 



V,^ = -^W2W^ + (3iqn 



= aw'' + ^g'^ 



(165) 
(166) 



E. A linearised model 



The final equations for the coupled three-component model are obviously rather complex. 
This is not surprising since, apart from working in a specific observer frame, we did not make 
any simplifications. Hence, the model is quite general, including the relevant nonlinearities 
and redshift factors. The main take home message should be that the steps involved in the 
derivation are natural and intuitive, but the expressions involved will be messy. One may 
query the immediate usefulness of the analysis, as it takes us far beyond what is currently 
considered in applications. However, the argument on behalf of the defence is clear. Once 
we have worked our way through the general case it is relatively straightforward to reduce 
the complexity by considering specific models. This is, in fact, a valuable exercise as it 
provides a clearer insight into the key features of the hot system. 

A natural, and in many cases of interest reasonable, assumption is that we only need to 
retain the linear relative velocities. As in the cold model, we neglect higher order terms in w°' 
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and q"", and we also ignore all the redshift factors by taking 7x ~ 1. It then follows naturally 
that the pressure is \E' = P and the temperature is = T. In the interest of clarity we will 
also ignore the resistive scattering between entropy (phonons) and the material components, 
i.e. we set Sy^ = 0. These assumptions lead to, i) the heat equation; 

K/3iqa+ (l + qa = V^T + Tii^ + l^2Wa + + + I32w') V.Mfe] , (167) 



ii) the generalised version of Ohm's law; 
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and, iii) the total momentum conservation equation; 



(168) 



(P + p) 11'+ X"^" VaP + g^Vau" + q'Vau''+ ±^ g'^ 
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= e (np - ne) i?" + e-^ (pp + Pe) e^-^^w^.P, . (169) 
P + p 

We simplify these relations further by noting that we can reinstate the assumption of 
charge neutrality, as the issues alluded to after eq. fll3ip originate from the redshift factors. 
Thus we let np = He, which leads to a number of simplifications (the arguments are the same 
as in the cold case). Focussing on the proton-electron plasma, we also assume that pe ^ Pp. 
It follows that the charge current is (again) given by 



(170) 



We can use this to write (11670 in the elegant form 



qa = -K[ ±1 VbT + Tila + 2u'VlbQa] 



(171) 



where 



Qa = Piqa + —Ja 



The thermal relaxation is encoded in this quantity. 



(172) 
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Turning to the momentum equation (11690 . we have 



(P + p) VaP + q'^VaU' + q'VaU'' + 

Finally, we find that Ohm's law simplifies to [t^ 



hac 



JaBr 



(173) 
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At this point, we have stripped the hot plasma model down to the level where it is easy 
to compare the final expressions to those of the cold model. At the linear level, the only 
difference is in the presence of the couplings that arise due to the entropy entrainment 
(expressed in terms of the different f3 coefficients) and the explicit presence of the heat 
flux in the momentum equation. These may seem like minor adjustments, but they are 
signiflcant. In particular, we need to retain the relevant relaxation times in order to ensure 
that the model is causal. 

Of course, the main differences between the two models we have developed enters at 
the nonlinear level (in the relative velocities). At quadratic order, the problem is non- 
adiabatic (as Fg 7^ 0) and it is no longer natural to assume charge neutrality. Given these 
effects, it would be very interesting to study a quadratic model in more detail. However, 
the corresponding problem is somewhat involved so we prefer to postpone discussion of it 
for the future. 

The model is completed by three scalar relations (whose origin are the conservation laws 
for the fluxes). In the linear case, these take the simple form 



p + {P + p)9 + Vaq^ - Ubq' = 



and 



which can be replaced by 



s + se + Vaij;] =0 , 



Var = 



he + rieO = . 



(175) 
(176) 

(177) 

(178) 
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V. CONCLUDING REMARKS 



We have developed the theory for charged fluids coupled to an electromagnetic field in 
the framework of general relativity, accounting for both a phenomenological resistivity and 
the relaxation times (associated with the charge current and the heat fiux) that are required 
to ensure causality. The final formalism can be applied to a range of interesting problems 
in astrophysics and cosmology. The cold two-component plasma model (from Section III) 
extends the "ideal" magnetohydrodynamics framework in several directions, and the hot 
model (from Section IV) adds dimensions that come into play when thermal aspects of the 
problem cannot be neglected. These developments are important as a number of interesting 
problems may require "non-ideal" aspects for their solution. Of most obvious relevance are 
problems involving not only electromagnetic fields but the live spacetime of general relativity. 
Several key gravitational- wave sources come to mind, like core-collapse supernovae ^| and 



compact binary mergers 
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47( 1 . Both cases involve strong gravity, a significant thermal 



component and magnetic fields. To apply a resistive framework to these problems is, of 
course, seriously challenging but this does not mean that we should not have aspirations in 



this direction 134 



48| are also of obvious relevance. 



36|. Actual multi- fluid simulations 
Focussing on relativistic stars, one can think of a number of unresolved problems, ranging 
from the dynamics of the magnetosphere and the pulsar emission mechanism to the formation 
and evolution of the star's interior magnetic field. These are problems where there has been 
significant progress, but further effort is required. In the case of the magnetosphere, the 
main focus has been on force-free models, but recent arguments 49| point to the need to 
include resistivity in the discussion. In the case of the formation and evolution of a compact 
star's global magnetic field, we need a better understanding of dynamo effects that may come 
,„to operation, see 0a.da.o [501 fo.. a. .ece,. review) and we also need to understand the 



51| . There are some 



coupled evolution of the star's spin, temperature and magnetic field 
very difficult issue to resolve here. 

In fact, the suggested examples highlight the need to develop the theory further. Typical 
questions that would need to be addressed involve (i) the dynamics predicted by the model, 
e.g. causality and stability of wave propagation and relation to issues hke pulsar emission 
or the launch of outfiows and jets, (ii) transitions between spatial regions where different 
simplifying assumptions are valid, such as a region in the magnetosphere where the fluid 
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model applies and a low density region where the description breaks down and one would 



need to fall back on a kinetic theory model 
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44| , the transition from magnetosphere to 



interior field at the star's surface or, indeed, accreting systems where an ion-electron plasma 
describes the inflowing matter while regions in the magnetosphere may still be appropriately 
modelled as a pair-plasma, (iii) the role of more complex physics, like the superconductor 



that is expected to be present in the star's core [52(] or regions where the assumption that 
the medium is electromagnetically "passive" does not apply, possibly in the pasta region 
near the crust-core transition. The present work provides a foundation for developments in 
all these directions, but each problem is associated with specific challenges that will need to 
be addressed if we want to make further progress. 
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npUp + rieV^ = . 
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readily extended to the "hot" case, where heat flux is added to the problem, as in Section IV. 
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for proton-electron plasmas. It is also clear that we have no reasons to assume that the 
gradients 

obey a similar ordering. 
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